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Abstract The Helioseismic and Magnetic Imager (HMI) aboard the Solar Dy¬ 
namics Observatory (SDO) observes the Sun at the Fe i 6173 A line and returns 
full-disk maps of line-of-sight (LOS) observables including the magnetic flux 
density, velocities, Fe I line width, line depth, and continuum intensity. These 
data are estimated through an algorithm (the MDI-like algorithm, hereafter), 
which combines observables obtained at six wavelength positions within the Fe 
I 6173 A line. To properly interpret such data it is important to understand any 
effects of the instrument and the pipeline that generates these data products. 
We tested the accuracy of the line width, line depth, and continuum intensity 
returned by the MDI-like algorithm using various one-dimensional (ID) atmo¬ 
sphere models. It was found that HMI estimates of these quantities are highly 
dependent on the shape of the line, therefore on the LOS angle and the magnetic 
flux density associated with the model, and less to line shifts with respect to the 
central positions of the instrument transmission profiles. In general, the relative 
difference between synthesized values and HMI estimates increases toward the 
limb and with the increase of the field; the MDI-like algorithm seems to fail 
in regions with fields larger than approximately 2000 G. Instrumental effects 
were investigated by the analysis of HMI data obtained at daily intervals for a 
span of three years at disk center in the quiet Sun and hourly intervals for a 
span of 200 hours. The analysis revealed periodicities induced by the variation 
of the orbital velocity of the observatory with respect to the Sun, and long-term 
trends attributed to instrument adjustments, re-calibrations and instrumental 
degradation. 
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1. Introduction 


The Helioseismic and Magnetic Imager (HMI) aboard the Solar Dynamics Ob¬ 
servatory (SDO) samples the magnetically sensitive photospheric Fe i 6173.34 A 
absorption line, through six narrow-band filters in right- and left- circular polar¬ 
ized light. Like the data-reduction pipeline of the Michelson Doppler Interfer¬ 
ometer (MDI) on board the Solar and Heliospheric Observatory (SOHO) (which 
sampled the Ni 6768 A line at five wavelength positions), the HMI-pipeline 
employs an algorithm (the MDI-like algorithm) that combines the filtergrams 
to obtain measurements of line-of-sight (LOS, hereafter) Doppler velocity, mag¬ 
netic flux density, and Fe i continuum intensity. As a by-product, the MDI-like 
algorithm also provides estimates of the Fe i line depth and width. We refer to 
the continuum intensity, line depth and width as line parameters, hereafter. 

The MDI-like algorithm assumes that the profile of the observed line is Gaus¬ 
sian shaped. Such an assumption, together with uncertainties in the transmission 
profiles of the instrument, saturation of the line in the presence of strong mag¬ 
netic fields, and Doppler shifts induced by plasma motion and solar rotation, 
are known to have caused uncertainties in measurement s retur n ed by the MDI. 
For in stance. Wachter, Schou. and Sankarasubramanian ( 20061 1 . Raiaguru et al. 
( 2007 1 and Wachterl (I2008ll investigated uncertaintie s in D opple r velocity es- 


timates. Tran et al. 


Demidov and Balthasar ( 2009l l and Ulrich et al. 


( 2 OO 9 II investigated uncertainties in ma gnetic flux measurements, while Mathew 


et al. (l2007l l and ICrisciioli et all (|201li l focused on uncertainties in continuum 
intensities. 

Similarly, recent studies have inve stigated the accuracy of HMI measurements. 
Fleck. Couvidat. and StrausI ( 20Illl employed three-dimensional hydrodynamic 
simulations to study the effects of uncertainties in the shap e and pos i tion o f 
the HMI transmission profil es on Doppler veloc ity estimates. Lin et al\ 1 2012h . 
Pietarila et all 1 2013h and Rilev et al\ 1 2014l l compared HMI magnetic flux 
measurements with those obtained with MDI and other instruments and derived 


conve rsion factors between the various magnetograms. The analysis bv iLiu et al 


(I 2 OI 2 I I showed in particular that HMI magnetic flux density measurements are 
affected by 24 and 12-ho ur periodicities, induced by the SDO orbital motion. 

Couvidat et al. ( 2012l l employed high spatial resolution spectro-polarimetric 
observations of an active region at disk center obtained at the Dunn Solar 
Telescope with the Interfe rometric Bidimensional Spectropolarimeter (IBIS) in¬ 
strument I Cavallini . 200^ to investigate the effects of line shape and velocity 
on all the data products (LOS velocity, magnetic flux density and Fe i line 
parameters) returned by the MDI-like algorithm. These authors found that the 
differences between measurements derived from the IBIS observations and those 
derived applying the MDI-like algorithm to the observed spectral profiles are in 
general below 20%, with the best estimates obtained for the continuum intensity 
(below 2%). 

Besides uncertainties resulting from assumptions in the derivation of the MDI- 
like algorithm, the HMI measurements are known to be affected by instrumental 
effects. These include a long-term increase of the opacity of the entrance window, 
variations of the focus, and uncertainties in the shape and position of the filter 
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transmission profiles (Couvidat (2014), private communication). These effects 
are partially compensated for in the case of velocity and magnetic field measure¬ 
ments along with continuum intensity measurements, but not so for the Fe i line 
depth and width. Moreover, the effects of variations of instrument characteristics 
on the HMI data product, especially on long-temporal scales (months to years) 
have not been presented in the literature yeta 

It is important to stress that the main purpose of the HMI-pipeline is to 
provide Dopplergrams, magnetograms, and Fe i continuum intensity, while the Fe 
I line width and depth are a by-product of the MDI-like algorithm. For this reason 
little effort has been dedicated to the improvement of such measurements (see 
also Section[2|). On the other hand, investigations of properties of line profiles are 
interesting for several types of studies, especially those in t he framework of sun- 
as-a-star and long-term variations of solar magnetism (e.q., Crisciioli et al . , 201 ,11 : 


Bertello. Pevtsov, and Pietarilal. 201^ PietariiaandljivlngstonL 2011 ^Livingston 

et ai, 200?!: Penza. Pietropaolo. and Livingstonl. ioOfih. 

In this study we therefore investigate the accuracy of HMI data products, 
focusing on the Fe i line-shape parameter estimates. To this end we employ 
synthetic Fe i line pro files to test the MDI-like algorithm, extending the work by 
ICouvidat et al. (20l3) to a larger sample of different LOS values, Doppler shift 
and magnetic field strength. We also employ HMI data to investigate the effects 
of SDO’s orbital motion and instrumental degradation on short (days) and long 
(months to years) term measurements of the Fe i line parameters. 

The article is organized as follows: the implementation of the MDI-like algo¬ 
rithm is briefly described in Section [21 tests of this algorithm using synthesized 
line profiles are discussed in Section [31 the analysis of HMI data is presented in 
Section |4l and our findings and their implications are summarized in Section |5l 


2. The MDI-like Algorithm 


The HMI uses a narrow band tunable filter to sample the Fe i 6173 A line at six 
different wavelength positions. 

Figure [T] displays an example of the transmission profiles for each of the 
wavelength positions alon g with a synthesized ir on line using a ID quiet-sun at¬ 
mosphere model (FAL-C) ( Fontenla et al llHii). As the line is sampled in both, 
left and right-circular polarized light (LCP and RCP, hereafter), 12 filtergrams 
are obtained during the wavelength tuning process. After standard data correc¬ 
tion (flat field, dark current, cosmic rays), these filtergrams are subsequently 
used to compute derived data products {i.e. observables) like Dopplergrams, 
magnetograms, line-width, line-depth and continuum-intensity maps employing 
a MDI-like algorithm. In synthesis, the algorithm estimates the data products 
as functions of the first two discrete Fourier coefficients of the Fe i line profile, 
which is assumed to be Gaussian-shaped. The Fourier coefficients are estimated 
by a proper combination of the filtergram intensities. A complete description of 


^Some information on instrument characteristics/calibrations can be found at 
http: //jsoc.stanford.edu/ 
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Figure 1. Example of the six HMI transmission filter profiles. The thicker line is a synthetic 
Fe I 6173 A line computed using a FAL-C (quiet Sun) model. 


the algorithm and of its implementation is given in Couvidat et al. ( 20121 '). In 
the Appendix we summarize only the steps necessary for the estimation of the 
Fe I line parameters. 

We implemented the algorithm in the Interactive Data Language (IDL). The 
main difference between our implementation and the one described in Couvidat 


et al. (120121) is in the HMI map employed to estimate the tr parameters used to 
compute the continuum intensity and line depth, as explained in the Appendix. 

Estimates derived by the MDI-like algorithm are affected by uncertainties 
induced by the physical and numerical assumptions on which the algorithm 
was derived. Namely, these uncertainties are generated by deviations of the Fe 
I line from a Gaussian profile, the deviation of the transmission profiles from 
A-functions, the assumption that a discrete Fourier transform can be computed 
based on only six data points, and the use of only two orders in the Fourier 
expansion. Moreover, changes of the shape of the transmission filter prohles, 
together with relative shifts of the transmission filter positions with respect 
to the center of the line (induced either by drifts of the transmission profiles 
or Doppler effects), also affects the estimates of the data products. The HMI- 
pipeline employs proper look-up tables to partially compensate for such effects. 
Nevertheless, these corrections are applied only to the Doppler velocity and mag¬ 
netic field and thus partially to the Fe i continuum intensity, which is function 
of the Doppler shift (see Equation (jH]) in the Appendix). The Fe i line depth 
and width estimates are a by-product of the MDI-like algorithm and no further 
correction is applied to them. 
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Table 1. ID atmosphere models used 
to synthesize the Fe i line and cor- 
resp onding vertical magnet ic field val¬ 
ues dFontenla et aZl Il99flri . The Ku- 
rucz models are denoted Kur_XXXX 
where XXXX is effective temperature in 
Kelvin. 


Model 

Description 

B [kG] 

FAL-A 

quiet Sun 

— 

FAL-C 

quiet Sun 

— 

FAL-E 

faint network 

— 

FAL-P 

facula 

1.00 

Kur_5500 

sunspot 

1.10 

Kur_5250 

sunspot 

1.70 

Kur_4750 

sunspot 

2.25 

Kur_4500 

sunspot 

2.50 

Kur_4250 

sunspot 

3.00 

Kur_4000 

sunspot 

3.50 

Kur.3750 

sunspot 

4.00 


3. Tests on Synthetic Profiles 


In order to investigate the effects of the deviation of the line shape from a 
Gaussian and the effects caused by the shift of the line with respect to the HMI 
wavelength filter positions we used synthesized profiles of the Fe i 6173 A line 
as an input to the MDI-like algorithm. The synthesized input profil es were com¬ 
puted in non-local thermodynamic equilibrium with the RH code (jUitenbroekl . 
200111 on a 3 A wide spectral range centered around the line. 

Since the shape of a spectral line is determined by the physical properties 
of the plasma, and, in the case of magnetically sensitive lines, by the magnetic 
field strength we synthesized the line using various static one-dimensional atmo¬ 
sphere models representing different features observed on the solar disk. Quiet 
sun, network and facu lar region profiles were synthesized using FAL99 models 
( Fontenla et al. . 1999l l. while sunspot profiles were synthesized using Kurucz 
models characterized by various effective te mperatures. We also p erformed tests 
on the three Malt by models as described in Criscuoli et all ( 2011 1. nevertheless 
results obtained with these latter models are qualitatively in agreement with 
those found from the Kurucz models and therefore they are not discussed here. 
To investigate the effects of the Zeeman splitting on lines, the facular model 
FAL-P and the sunspot models were synthesized imposing a vertical magnetic 
field of constant strength. Although the imposed magnetic fields are arbitrarily 
associated with the model effective temperature, the values do follow the widely 
known anti-c orrelation between magnetic field and temperature in sunspots {e.g., 
SchadL 20141 . and references therein). The models employed and the correspond¬ 
ing magnetic field strengths are summarized in Table [1] For each model we 
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Figure 2. Examples of Fe I line profiles synthesized through three ID atmosphere models at 
three viewing angles. Continuous line: FAL-C. Dashed line: FAL-P with uniform magnetic field 
of 1000 G. Dot-dashed line: Kurucz_4250 with uniform magnetic field of 3000 G. All profiles 
have been normalized to the corresponding continuum intensity at disk center. 





aa[A] aa[a] aa[A] 


considered eight viewing angles, expressed in the following as the cosine of the 
heliocentric angle, /i. Note that by synthesizing the line for different viewing 
angles we equivalently investigate the effect of inclined magnetic fields. Examples 
of line profiles generated using these models are shown in Figure[2] To investigate 
the effects of velocities due to convective motions, solar rotation and relative 
motions between the observatory and the Sun, we also shifted the line profiles 
between [-4,4] km s“^ (note that the SDO velocity relative to the Sun spans 
between -3.2, and 3.2 km s“^). The line parameters were estimated by feeding 
the synthetic profiles to the MDI-like algorithm and the obtained values were 
compared with those derived from direct estimates on the synthetic profiles 
(synthetic values, hereafter). In particular, following the definitions of the MDI- 
like algorithm (see Appendix) for the line parameters, we defined the synthetic 
continuum as the intensity synthesized at 0.5 A from line center, the line depth 
as the difference between the emergent intensities in the continuum and at the 
center of the line, while the line width was estimated as the full width at half 
maximum. All line parameters were derived separately for the TCP and RCP 
components of the spectral line and then averaged to give the final value. 

3.1. Line Shape Effects: Model and Line-of-Sight Dependence 

Figure [3] shows the relative difference amongst the line parameters derived with 
the MDI-like algorithm and the synthetic ones, for the FAL99 models and for 
various lines-of-sight. The plots show that the relative difference obtained from 
FAL-A, FAL-C, and FAL-E are similar and that in particular the MDI-like 
algorithm on these three models overestiiuates the line width by 5-10% and 
the line depth by 10-30%, with a higher relative difference in the line depth 
resulting at shallower lines-of-sight (smaller fi values). Continuum intensity is 
found to be more accurate as the difference between synthetic values and MDI- 
like algorithm estimates are within 1% with a minimum of about 0.1% close to 
disk center. The differences obtained with the FAL-P model representative of 
a facular region with a vertical magnetic field of 1000 G, are larger and have 
opposite sign with respect to the other FAL models investigated. This is because 
the Fe i line shape profiles obtained from the first three models are similar and 
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Figure 3. Relative differences between synthetic and MDI-like algorithm-calculated line pa¬ 
rameter values for the FAL99 models: FAL-A (open circle), FAL-C (open square), FAL-E 
(X), and FAL-P (open triangle). FWHM denotes the line-width, LD the line-depth and IC 
the continuum intensity. The bar on top of the acronyms denotes quantities returned by the 
MDI-like algorithm. 


close to the shape of a Gaussian profile whereas the shape obtained from the 
FAL-P is broadened by the Zeeman effect (as shown in Figure^). The line width 
is in fact underestimated by 15-20%, the line depth is overestimated by 20-60%, 
and the continuum intensity is underestimated by 0.5-2%. The differences for 
this model decrease toward the limb, because, as shown in Figured the line core 
intensity enhancement due to the Zeeman splitting decreases for shallower LOS. 

The results obtained for all seven Kurucz models are plotted in Figured) The 
relative difference between the synthetic values and those returned by the MDI- 
like algorithm are very large, with the lowest values found in the case of the 
continuum intensity, for which the discrepancy is up to 20%. For the line depth 
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Figure 4. Relative differences between theoretical and MDI-like algorithm-calculated line 
shape values for the Kurucz sunspot models: Kur_5500 (open circle), Kur_5250 (open upward 
triangle), Kur_4750 (open diamond), Kur_4500 (open square), Kur_4250 (star), Kur_4000 (X), 
Kur_3750 (open downward triangle). FWHM denotes the line-width, LD the line-depth and 
IC the continuum intensity. The bar on top of the acronyms denotes quantities returned by 
the MDI-like algorithm. 


and the line width the differences are up to 100%. In general, the differences in¬ 
crease with the increase of the magnetic field strength, as a consequence of both, 
the line saturation and the poor sampling of the wings and continuum in the case 
of large line broadening. In particular, for magnetic field intensities larger than 
approximately 2000 G the MDI-like algorithm returned unphysical line width 
estimates. These cases are marked above the double line in the plots. This occurs 
because for these models the sum of the squares of the second Fourier coefficients 
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is greater than the sum of the squares of the first Fourier coefficients, so that the 
argument of the logarithm in Equation Q in the Appendix returns a negative 
number whose square root is not real. Similarly, the ’infinite’ relative difference 
found at vertical lines-of-sight for the line depth (these cases are also beyond the 
double line in the plot) are due to the saturation of the line at large magnetic field 
strengths (see Figured]), so that the synthetic line depth tends to zero, while 
the MDFlike algorithm, expecting a Gaussian-shaped profile, returns a larger 
value. It is also interesting to note that for the largest magnetic field strengths 
investigated the deviations from the synthetic values of the line depth and of the 
line width present small center-to-limb variation because of saturation effects. 
For the same reason, as already pointed out in previous studies investigatin g the 
accuracy of the MDI algorithm ( Mathew et al. , 2007 ; Criscuoli et all 201 ill , the 
discrepancies of the continuum intensities do not increase monotonically with 
the increase of the magnetic field strength, instead they decrease for the largest 
field values investigated. 


3.2. Line-Shift Effects 


We also investigate the influence of a line shift on the determination of the line 
parameters. In the following we report only results obtained with the FAL99 
models, as those obtained for sunspots are qualitatively in agreement. Figure 
[S] shows the relative difference between the values obtained from the MDI-like 
algorithm with the line at rest and the line shifted versus the amount of shift. 
The plots reveal a periodicity for all the line parameters, which is propo r tional 
to the spacing betwe en the HMI transmission profiles (see also IWachteii . l2008l : 


Criscuoli et al . 201l[ l.. 


The variations obtained for models FAL-A, -C, and -E are similar in amplitude 
(a few percent) and sign, whereas the results obtained with the FAL-P model 
show larger deviations (up to 13%) and opposite sign, in agreement with the 
results presented in Figure [3| The increase of the difference between measured 
and synthetic values is due in this case to the fact that a line shift causes the left 
and right wings of the line to be unevenly sampled. The broadening induced by 
the Zeeman splitting enhances this effect (for larger broadenings even modest 
shifts may cause portions of one of the wings not to be sampled at all), so that 
it is no surprise that the uncertainties increase with the increase of the magnetic 
held strength associated with the model. 


In agreement with results presented in ICriscuoli et al\ (|201lh for the contin¬ 
uum intensity estima ted by the MDI algorit hm, the plots also show asymmetries 
around zero velocity. Criscuoli et all ( 201lll ascribed this to the fact that a line 
shift could compensate for or enhance the original shift of the lines that they 
employed (obtained from non-static simulations and observations) for their tests. 
Because the results presented in Figure[S]have been obtained with static models, 
the asymmetries around the zero velocity must be ascribed to small asymmetries 
in the hlter transmission prohles. 

Finally, it is worth to note that the results reported in Figure [5] allow to 
address effects introduced by line shifts only. In practice, uncertainties are in¬ 
troduced by the combined effect of line-shape and line-shift. For each model and 
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Figure 5. Relative differences of line parameters calculated with the MDI-like algorithm 
between Doppler shifted line and line at rest, for the Fe i line generated by the Fontenla models 
at fj. = 1. Solid line: FAL-A. Dashed line: FAL-C. Dot—dashed line: FAL-E. Triple-dot—dashed 
line: FAL-P. 


shift, it is easy to show that, if h is the uncertainty resulting from the deviation 
of the line shape from a Gaussian profile (results reported in Figure [3]), and h' 
is the uncertainty resulting from a line-shift (results reported in Figure [5]), then 
the uncertainty resulting by the combination of the two effects is h + h' + h ■ h'. 


4. Tests on HMI Data 


The results presented in the previous section indicate that the accuracy with 
which the MDFlike algorithm estimates the Fe i line parameters strongly de¬ 
pends on the shape of the line profiles and the line shift. The first effect domi¬ 
nates, at least for shifts in the range investigated. We expect studies aiming at 
investigating the center-to-limb variation or th e tempora l evolu tion of magnetic 
features to be particularly affected. For instance lLiu et al. (2012) showed that the 
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evolution of the LOS magnetic flux density of an active region (AR) measured 
by HMI presents 12 and 24 hours cyclic variations, which are generated by the 
orbital motion of the SDO. This signal is super-imposed to a longer temporal 
trend resulting by the combination of the change of the LOS with solar rotation 
Doppler effects as the region crosses the solar disk. To investigate these effects 
also on the lin e-shape paramet ers returned by the HMI, we repeated the analysis 
performed bv iLin et al. 1 2012h on a different region (i.e. AR 11092) and inves¬ 
tigated the evolution of the line parameters and of the magnetic flux density on 
quiet, facular, and sunspot regions. 

We also analyzed the variation of the line-shape parameters over a data set 
of HMI observations spanning about three years of observations with the aim 
to investigate the effects of instrumental degradation described in Section [TJ As 
properties of magnetic features can vary over the cycle, for this last analysis only 
quiet-sun regions at disk center were considered. 


4.1. Short-Term Line Parameters Variations 


In order to investigate short-term effects of solar rotation and orbital velocities 
on HMI data products we analyzed a one hour cadence data set spanning about 
200 hours in total, acquired between July 30 and August 8, 2010. In particular, 
we analyzed the 720 sec magnetic flux, line-depth, continuum intensity and 
line-width data products. We employed magnetograms corrected for the fore¬ 
shortening {B = Bobs/fJ-) to identify quiet (pixels where absolute magnetic flux 
density is < 100 G), facular (pixels where absolute magnetic flux density is > 
100 and < 1000 G) and sunspot (pixels where absolute magnetic flux density 
is > 1000 G) regions around AR 11092 (NOAA number). This active region 
was chosen because it showed little evolution during its passage over the disk. 
For each class of pixels we then studied the temporal variation of the average 
Fe I line parameters and the magnetic flux density. The investigation of the 
average magnetic flux density allowed verifying that the observed line param¬ 
eter variations were not caused by the evolution of the magnetic field in the 
region. The variations of the average magnetic flux density computed over the 
three classes of regions are reported in Figure [5] (left panel). The dashed and 
continuous lines represent the magnetic flux density non-corrected and corrected 
for the for e -short ening, respectively. These plots agree with results reported in 
Liu et oZI ( 2012f l. which showed that the magnetic flux density of an active 
region non-corrected for the fore-shortening is maximum at disk center, while 
the flux density of quiet regions does not depend on the position on the disk. 
When correcting for the fore-shortening, the curves change concavity, with a 
minimum^obseiwed_a4_dis^,centerj,^ug£estin^^ha^^he,^rends,j;e2orted_j^Liu 
et al. (|2012h are mostly caused by projection effects. Semi-diurnal a nd diurnal 
period ic variations of the flux are also observed. These, as explained in lLiu et all 
(1291211 . can be mostly ascribed to spacecraft orbital motions effects not entirely 
compensated for by the look-up tables. The mostly symmetric trends around 
disk center suggest that the magnetic field of the region stayed fairly constant 
during the observation time. 

The right panel in FigurelHlshows results obtained for the line parameters. The 
plots show that, like the magnetic flux density, line parameters are symmetric 


SOLA: hini_article_v3.tex; 10 Februciry 2015; 1:40; p. 11 






















sunspot [G] B-foculae[G] B-quiet [G] 


Cohen et al. 






Figure 6. Variation of average magnetic flux density (left) and average line shape parameters 
(right) of quiet, facular and sunspot pixels (see text) observed around AR11092 by the HMI. 
The observations started on July 30, 2010 at 11:58:26 UT. The region was at the closest 
proximity to disk center about 100 hours after the beginning of the observations. Continuum 
intensity and line depth are normalized to the average continuum intensity of quiet pixels at 
disk center. 


around the time of closest proximity to disk center. These variations must be 
mostly ascribed to the combination of LOS effects and solar rotation. The typical 
periodicities caused by orbital motions of the satellite are also observed. These 
are more noticeable in the case of the line width in magnetized pixels, as expected 
by results discussed in Section It is important to notice that all the observed 
temporal variations are due to both LOS and velocity effects and that it is 
not possible, in practice, to distinguish between them. However, comparison of 
results presented in Section iTTl and Section|32]suggests that the longer temporal 
trends of the continuum intensity and line depth reported in FigureElare mostly 
caused by LOS effects rather than by the solar rotation, because variations 
induced by this last effect (the solar rotation is ± 2 km s“^ at the equator) are 
about one order of magnitude smaller than variations induced by LOS effects. 
As a consequence, because, as shown in Section IXTl the shape of a line in high 
magnetized regions is dominated by the Zeeman splitting rather than by LOS 
effects, the line parameters show almost no variation in sunspot pixels. 

4.2. Long-Term Line Parameters Variations 

In order to investigate long-term variations of the Fe i line parameters returned 
by the HMI, we have analyzed a set of magnetic flux density, line-depth, con- 
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tinuum intensity and line-width data acquired daily between May 1, 2010 and 
May 30, 2013. For this analysis we also employed 720 sec data, obtained as 
closest as possible to 12:00 UT. In order to more easily distinguish between 
effects inherent to the measurements from those induced by variations of mag¬ 
netic feature properties induced by the solar magnetic cycle, we restricted the 
analysis to quiet pixels located close to disk center. To this end we analyzed 
the temporal variation of the average value of the Fe i data products computed 
over pixels where the absolute magnetic flux density is < 100 G, located in 
a square area of 500x500 pixels (about 250x250 arcsec) around disk center. 
The results are shown in Figure [T] The line-shape parameter values have been 
normalized to the value measured the Hrst day analyzed. The plots of line-shape 
parameters show an overall decrease with time and several discontinuities. The 
long-term trends are caused by various issues, mentioned in Section [1] that 
are known to affect the HMI, such as an increase of the opacity of the entrance 
window (most likely induced by exposure to UV radiation), changes in the shape 
and positions of the transmission profiles and variations of the focus caused by 
temperature changes in the entrance window especially after eclipse periods. 
The discontinuities are due to periodic re-calibration of the instrument that 
are meant to mitigate these effects. The re-calibrations include: an increase of 
the standard exposure time, retuning of the transmission profiles (the days of 
the retuning are marked with arrows in the figure), adjustments of the focus, 
and front window re-heating after eclipse periods (calibrations calendar and 
descriptions are available at http://jsoc.stanford.edu/doc/data/hmi/). Among 
these effects, those that dominate the long-term trends are the increase of the 
entrance window opacity and the shift of the transmission profiles. Comparison 
of plots in Figure [7] and Figure [^indicate that the continuum intensity and line 
depth are mostly affected by the increase of the entrance window opacity rather 
than by a shift of the transmission profiles. The observed long-term decrease of 
these quantities is in fact about 10%. Plots in Figure [5] show that in the case of 
the line depth such variation can be generated by shifts larger than about 4 km 
s“^ (which seems an unlikely large value, considering that, as stated above, the 
transmission profiles are periodically re-tuned), while the continuum intensity is 
almost insensitive to relative shifts of the line with respect to the positions of 
the transmission profiles. Instead, Equation Q in the Appendix indicates that 
the line-width is insensitive to intensity variations, as long as these variations 
are the same for all the six filtergrams. The observed decrease of about 4% for 
this parameter is compatible with shifts of the transmission profiles smaller than 
1 km s“^. 

It is also noteworthy that the observed variations of the line parameters cannot 
be ascribed to variations of the quiet-Sun magnetic field. In fact. Figure [7] shows 
that the average absolute magnetic flux density slightly increases with time, 
which would cause both the line width and depth to increase. Note also that, 
as mentioned above, the magnetic flux density measurements are compensated 
for line shift effects (induced by both LOS orbital velocity changes and drifts 
of the transmission profiles) and that the magne tic flux density i s inde pendent 
from the absolute intensity of the filtergrams (see ICouvidat et al . , 2012 h , so that 
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Figure 7. Line parameters and average absolute magnetic flux density observed by the HMI 
in quiet pixels close to disk center over a period of three years. The arrows indicate the days 
in which the HMI filters were re-tuned. 


the small flux-density increase over time, in agreement with the increase of the 
magnetic activity during the period analyzed, could be real. 

Finally, the plots in Figure [7] show periodic variations of the line parameters, 
which are particularly noticeable in the case of the line width. A comparison 
of the Fourier spectra of the SDO velocities and of the line width (not shown) 
confirmed that even these oscillations are due to the relative variation of the 
SDO orbital velocity with respect to the Sun. 


5. Discussion and Conclusion 

In this study we investigated the accuracy of the Fe i 6173 A line-parameter 
values {i.e. line depth, line width, continuum intensity) returned by the HMI. 
To this aim, we first tested the accuracy of estimates returned by the MDFlike 
algorithm, which is part of the HMI data reduction pipeline. We then identified 
systematic and instrumental effects by studying the temporal evolution of HMI 
data obtained at different cadence over short (days) and long (years) temporal 
ranges. 

An MDI-like algorithm was implemented and tested with synthetic Fe i 6173 
A line profiles generated using ID atmosphere models representing quiet-sun, 
faint network, facular, and sunspot regions at different viewing angles across the 
solar disk. We found that the accuracy of this algorithm strongly depends on 
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the line shape determined by magnetic field strength and viewing angle. The 
algorithm works relatively well in quiet and facular regions, where the accuracy 
of the continuum intensity is better than 2% and the accuracy of the line depth 
and line width are on average of the order or below 20%. The accuracy decreases 
with the increase of the magnetic field strength, so that in highly magnetized 
regions (e.g. sunspots) it is in general of the order of several tens of per cent, 
with the best accuracy found again for the continuum intensity (about or below 
20%). For a magnetic field strength larger than about 2000 G the algorithm 
encounters numerical problems due to the saturation of the line, so that in these 
regions HMI results cannot be c onsidered reliab l e. Th ese results are in partial 
agreement with those obtained by Convidat et al. ( 2012h . who employed observed 
spectra to investigated the accuracy of the data products returned by the MDI- 
like algorithm. In particular, the accuracies that we estimated for the continuum 
intensity along vertical lines-of-sight are in good agreement with those presented 
by those authors. In the case of line width and depth, the agreement is still 
good in quiet and facular regions, but in regions of stronger magnetic fie ld we 
obtained a smaller accuracy than the one reported bv ICouvidat et ali (|2012l l. The 
discrepancy between our results and those presented by those authors has to be 
ascribed in part to uncertainties in the radiative synthesis and assumptions in our 
models {e.g. the atomic transitions values employed, the arbitrary association 
of an atmosphere model with a magnetic field strength, the assumption of a 
vertical constant field, and the non-inclu sion of molecules in umbra models). On 


the other hand, IConvidat et al\ (l2ni2n estimated a better accuracy in highly 


magnetized regions because they employed observed spectra on which the line 
shape variations induced by the Zeeman splitting (especially the increase of 
the core intensity) are strongly reduced in amplitude by scattered light (both 
sp atial and spectral) and finite spectral resolution effects (see also the discussion 
in Criscnoli et '^. l2mih . 


The effects of the line shift with respect to the nominal wavelength position 
of the filter transmission profiles have also been investigated, and they resulted 
to be in general smaller than LOS and magnetic field effects, at least for shifts 
in the range [-4,4] km s“^. 

Along with testing the HMI pipeline, we studied Fe i line width, line depth, 
and continuum intensity HMI data generated at both hourly intervals over a 
span of 200 hours in an active region and at daily intervals for a span of three 
years in quiet regions at disk center. We found both short-term and long-term 
periodicities caused by the orbit of the SDO, long-term decrease in intensity 
measurements due to opacity of the front window increasing, long-term decrease 
in line width due to drift in wavelength of transmission filters, along with 
discontinuities present in the long-term data due to instrument adjustments 
and re-calibrations. The LOS magnetic field and the LOS Doppler velocity are 
partially corrected for such effect, while the Fe i line parameters data are not. 
This complicates the analysis of temporal variations of the iron line parameters. 

Thus, the Fe i line parameter estimates returned by the HMI are affected 
by uncertainties due to approximations on which the MDI-like algorithm works 
(mainly, the assumption that the line is Gaussian shaped) and to instrumental 
effects. These uncertainties could be reduced with the use of proper look-up 
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tables, like those employed to correct LOS magnetic field and velocity estimates 
and thus continuum intensity measurements. However, such multi-dimensional 
look-up tables would likely be quite difficult to implement. Alternatively the 
algorithm might benefit from bei ng modified as to fit th e actual Fe i line profile 
rather than assume a Gaussian ( Couvidat et al . 2012h . In this case the speed 
of the fitting algorithm would be a decisive measure of its utility, as it would 
require fitting a very large number of pixels. The long-term effects induced by 
instrumental degradation, especially those induced by the increased opacity of 
the entrance window, are more difficult to correct for. It is possible to make 
such corrections if it is assumed that the quiet Sun does not vary, perhaps 
committing a small error, but long-term studies of the quiet Sun would therefore 
be comprised. Because such corrections are not available yet, results presented 
in this study are important for the correct interpretation of the Fe i 6173 A 
line-parameter estimates returned by the HMI. 
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Appendix 


A. The MDI-like algorithm 


We i mplemented the MDI-like algorithm following the description in Couvidat 
et al. (|2012i ). Here we summarize only the steps necessary to estimate the Fe i 
line parameters. 

The algorithm assumes a Gaussian shaped Fe i line profile: 


/(A) = r 


Id exp 
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rr2 


( 1 ) 


where Ic is the continuum intensity, Id is the line depth, and cr is a parameter 
that dehnes the width. Note that the algorithm returns as an estimate of the 
line width, the full width at half maximum (FWHM), which is proportional to 
(7 through the relation FWHM = 2yj2 In (2)tT. 

The algorithm first computes discrete approximations to the first and second 
Fourier coefficients using the six filtergram intensities Ij : 
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In our implementation the filtergram intensities were calculated as the convo¬ 
lution of the synthetic iron line with each of the six transmission filter profiles. 
The line amplitude and depth and the continuum intensity are then calculated 
according to the formulas: 


a = 



Id = 


T 


2a^/TT 


\lal+bl exp : 


(4) 


(5) 


1 

j=o 


Id exp 


(A,-Ao)2 


( 6 ) 


Here T = 412.8 mA is six times the nominal wavelength separation between the 
filter transmission profiles and Xj corresponds to the central wavelength of the 
j-th filter. 

The actual implementation of the MDI-like algorithm is slightly different than 
as explained above. Tests performed on Gaussian profiles have shown that the 
line parameters calculated using the formulas above are affected by inaccuracies 
generated by the finite term Fourier expansion (using only six wavelengths to 
calculate the Fourier transform), the finite sampling of th e line and the fa c t that 
the transmission profiles are not (5-functions. In particular, Couvidat et al\ ( 2012l l 
reported that Equation (|T|) overestimates the real width by ^ 20% and that 
Equation (O underestimates the actual depth by ~ 33%. To take into account 
these uncertainties, in the actual implementation Equation 0 is multiplied by 
a factor of 5/6 and Equation ([5|) is multiplied by a factor of 6/5. Tests have 
also revealed that the error in the continuum intensity is only ~ 1% so that no 
correction factor is applied to Equation (jH]). 

Moreover, Equations ([S]) and ([5]) are not computed using the u values returned 
by Equation ©• The reason is that, as shown in previous sections, the a returned 
by Equation (|4]) suffers from uncertainties induced by strong deviations of the 
line shape from a Gaussian profile (like in presence of strong fields), which would 
cause uncertainties in the estimate of the other two line parameters. Instead, 
the algorithm uses standard values of a derived from a HMI full-disk line-width 
map observed in a day of low activity. After the line width is converted to ct, the 
standard a values are computed as the fifth order polynomial fit to the variation 
of the azimuthally averaged (T-values with radial distance from disk center. 

The values of a as function of the cosine of the heliocentric angle employed in 
our implementation are illustrated in Figure [5] and were obtained from a HMI 
line-width map taken on September 10, 2009. 
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Figure 8. Line-of-sight dependence of the cr values employed to compute the line depth and 
the continuum intensity. The ic-axis is the viewing angle fi = cos 6 and the dashed line is a 
fifth-order polynomial fit to the data. 


SOLA: hini_article_v3.tex; 10 February 2015; 1:40; p. 19 




